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Background: The neutron skin of a heavy nucleus as well as many neutron-star properties are 
highly sensitive to the poorly constrained density dependence of the symmetry energy. Purpose: 
To provide for the first time meaningful theoretical errors and to assess the degree of correlation 
between the neutron-skin thickness of ^"^Pb and several neutron-star properties. Methods: A 
proper covariance analysis based on the predictions of an accurately-calibrated relativistic functional 
"FSUGold" is used to quantify theoretical errors and correlation coefficients. Results: We find 
correlation coefficients of nearly one (or minus one) between the neutron-skin thickness of ^"*Pb 
and a host of observables of relevance to the structure, dynamics, and composition of neutron stars. 

rsj Conclusions: We suggest that a follow-up PREX measurement, ideally with a 0.5% accuracy, 

,-H could significantly constrain the equation of state of neutron-star matter. 

o 

px| PACS numbers: 21.65. Cd, 21.65.Mn, 26.60.Kp, 26.60.-c 
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I. INTRODUCTION 



The equation of state (EOS) of neutron-rich matter plays a crucial role in understanding interesting phenomena in 

(— H both nuclear physics and astrophysics, such as the limits of nuclear existence, the dynamics of heavy ion collisions, the 

-^-^ structure of neutron stars, and the mechanism of core-collapse supernovae. A particularly important property of the 
1—1 EOS is the nuclear symmetry energy, whose value is well constrained by the ground-state properties of finite nuclei 

y only near saturation density [T] . Given that the symmetry energy reflects the increase in the energy of the system 
as protons are turned into neutrons (or viceversa), many nuclear and astrophysical observables are sensitive to its 
density dependence. For example, the EOS of neutron-rich matter is the sole feature responsible for the structure of 
neutron stars. Moreover, the neutron-skin thickness of heavy nuclei — a system that is 18 orders of magnitude smaller 
than a neutron star — is also sensitive to the symmetry energy. This is because the slope of the symmetry energy is 
related to the pressure exerted by the neutrons in creating both a neutron- rich skin and the neutron-star radius [H [3] . 
^^ Although considerable effort has been devoted to the understanding of the EOS of isospin asymmetric matter and 
(^ its possible impact on both laboratory and observational data [H [S] , our knowledge of the density dependence of 
^^ the symmetry energy remains incomplete. Whereas some of the laboratory /observational data is relatively easy to 

-y-! collect, most requires great perseverance and ingenuity. Thus, in an effort to identify observables sensitive to the 
f— ^ density dependence of the symmetry energy and to establish correlations among them, we rely on a powerful and 
^Sj systematic covariance analysis. Such a systematic statistical analysis may be used to attach meaningful uncertainty 

T— I estimates to theoretical predictions [B]. Moreover, it can be used to identify observables that, although at present may 

be beyond experimental/observational reach, display a strong sensitivity to the symmetry energy. Finally, through 
such a covariance analysis one can quantify the degree of correlation between various observables. It is the aim of 
the present contribution to establish for the first time quantitative correlations between the neutron-skin thickness of 



H 208pj-, g^j^j g^ variety of neutron-star observables. 

In a recent publication we developed a covariance analysis within a class of relativistic mean- field models [7] . Start- 
ing from a x^-minimization procedure, we presented a step-by-step implementation of the statistical approach that 
provided quantitative uncertainties in our theoretical predictions as well as robust correlations among physical ob- 
servables. In the present contribution we extend such a study to: (a) quantify the degree of correlation between the 
neutron-skin thickness of ^"^^Pb and several neutron-star properties and (b) provide meaningful theoretical uncertain- 
ties that arise from our incomplete knowledge of the density dependence of the symmetry energy. We have selected 
the neutron-skin thickness of ^''^Pb as it represent a laboratory observable of paramount importance in constraining 
the density dependence of the symmetry energy. Indeed, the strong correlation between the neutron-skin thickness 
of ^"^Pb and the slope of the symmetry energy at saturation density is a well established fact [H ISUll). Further, 



'Electronic address: 'Farrooh'Fattoyev@taniu-conimerce.edu| 
^Electronic address: jpiekarewicz@fsu.edu 



the highly anticipated Lead Radius Experiment (PREX) has just provided the first model-independent evidence of 
the existence of a significant neutron skin in ^"^Pb [12]. Building on the strength of the enormously successful 
parity-violating program at the Jefferson Laboratory, PREX used parity-violating electron scattering to determine 
the neutron-skin thickness of ^ospj^ ^q |-,g. 

-Rskin = ^ri — ^p = 0.33_Q;;^g fm , (1) 

where Rn{Rp) denotes the root-mean-square neutron(proton) radius. Although PREX achieved the systematic control 
required to perform this challenging experiment, unforeseen technical problems resulted in time losses that significantly 
compromised the statistical accuracy of the measurement. Thus, rather than achieving the original goal of a 1% error 
in the neutron radius, PREX had to settle for an error almost three times as large. Even so, the PREX measurement 
has already been used to impose some (mild) constraints on accurately-calibrated models of nuclear-structure [13j . 
Moreover, the — already approved — follow-up PREX measurement designed to achieve the original 1% goal will be of 
great value, especially if the unexpected large central value remains unchanged. 

Since the strong correlation between the neutron-skin thickness of ^^^Pb and the pressure of pure neutron matter 
near saturation density was first inferred [H H] , a large number of additional correlations between i?skin and several 
neutron-star properties have been uncovered [1|. These include: (a) the crust-to-core transition density [5], (b) 
neutron-star radii [H [T3], (c) threshold density at the onset of the direct Urea process [TS], and (d) the crustal 
moment of inertia [H [16], among others. Although most of these correlations have been previously established 
through a systematic variation of a previously calibrated model, in the present manuscript we attempt for the first 
time to reliably quantify all these correlations via a legitimate covariance analysis. A covariance analysis represents 
the least biased and most exhaustive tool to uncover correlations between physical observables [7[ [T7] . 

The manuscript has been organized as follows. In Sec. [TT[we briefly review the formalism required to implement the 



covariance analysis developed recently in Ref. [7^. In Sec. Ill we present physical arguments in favor of the expected 
correlations between the neutron-skin thickness of ^ospj^ g^j^^j several neutron-star properties. Properly estimated 
theoretical errors and correlation coefficients are also presented in this section. Finally, conclusions and suggestions 
for future work are summarized in Sec. IIVI 

II. FORMALISM 

The covariance analysis mentioned above will be implemented using the predictions from the accurately calibrated 
FSUGold model [TH] . The interacting Lagrangian density for the model — with its predictions generated at the rela- 
tivistic mean-field (RMF) level — is given by 



i?int = ^ 



ffs0- (.9vVp + |x • b^ + ^(l+r3)A^) 7^ 



v- 



This model contains an isodoublet nucleon field {ip) interacting via the exchange of two isoscalar mesons, the scalar 
sigma {(/)) and the vector omega {V'^), one isovector meson, the rho (b''), and the photon {A'^). In order to improve 
the standing of the model the Lagrangian density is supplemented by scalar and vector self-interactions. In particular, 
the scalar self-interaction (with coupling constants k and A) is responsible for reducing the compression modulus of 
nuclear matter from the unrealistically large value of K — 545 MeV obtained with the original Walecka model [111 [20] 
to about K — 230 MeV. Such a significantly lower value for the compression modulus is demanded by measurements 
of giant monopole resonances in medium to heavy nuclei |21H24[ . The quartic isoscalar- vector self- interaction is 
responsible for softening the EOS at high densities. Indeed, by tuning ^ one can generate different limiting neutron- 
star masses without modifying the behavior of the EOS around saturation density |25| . As such, C is fairly insensitive 
to laboratory observables and must be constrained from astrophysical observations. Finally, the mixed quartic vector 
interaction (as described by parameter Av) was introduced to modify the density dependence of symmetry energy [5], 
which is fairly stiff in the original Walecka model. By calibrating the parameters of the model using both ground-state 
properties of finite nuclei as well as their collective excitations [TS] , the FSUGold model has been fairly successful when 
compared against theoretical, experimental, and observational constraints |26j . However, as more observational and 
experimental data become available — such as the recently reported 2-solar mass neutron star [27] — further refinements 
may be required ,28j (for example by tuning Q. Yet, for the purpose of this contribution we will be satisfied with the 
use of FSUGold model as the basis for the covariance analysis. 

The covariance analysis follows closely the approach outlined in Refs. [T] [17] which, in turn, is based on the 
comprehensive text by Brandt [29j . The derivation, implementation, and power of the formalism was illustrated in 



our recent study using two relativistic mean-field models fT. Thus, in the present contribution we only offer a brief 
summary. The covariance analysis starts with the definition of the quality measure y^ . That is, 

where TV denotes the total number of selected observables that will be employed in the calibration procedure. Each of 
the observables. On , is assumed to have been determined experimentally with an accuracy of AO„, which acts as 
a weight factor in the quality measure. In addition, each of the iV observable is computed within the model On (p) 
as a function of the F model-parameters p= (pi, . . . ,p_f). The calibration procedure terminates when a set of optimal 
parameters po are found that minimize the quality measure. 

Given that our main goal is not the minimization procedure, but rather the estimation of theoretical uncertainties 
and the assessment of correlations among observables, we use a set of nuclear properties generated directly from 
the FSUGold model. This guarantees that the quality measure so defined is minimized. In particular, the following 
set of A^ = 8 observables have been selected as input for y^: (1) the saturation density poi (2) the binding energy 
per particle of symmetric nuclear matter at saturation density £o, (3) the binding energy per particle of symmetric 
nuclear matter at twice saturation density e(2po), (4) the compression modulus of symmetric nuclear matter Kq^ (5) 
the effective (Dirac) mass of symmetric nuclear matter at saturation density Mq , (6) the symmetry energy J evaluated 
at a sub-saturation density of pwO.lfm" , (7) the slope of the symmetry energy at saturation density L, and (8) 
the maximum neutron-star mass Afmax- We attach a 2% uncertainty to all of the observables — except for the case of 
the slope of the symmetry energy L where the significantly larger value of 20% is assumed. Such a large uncertainty 
reflects our poor understanding of the density dependence of the symmetry energy. Note that our choice of J, namely, 
the symmetry energy at sub-saturation density rather than at saturation density, follows from the fact that heavy 
nuclei constrain the symmetry energy at a density which is intermediate between the center and the surface of the 
nucleus. Indeed, at such a sub-saturation density the theoretical uncertainties appear to be minimized [1]. 

Having obtained the optimal parameter set po through the minimization of the quality measure, one can the proceed 
to compute and diagonalize the symmetric matrix of second derivatives. This matrix contains all the information 
about the behavior of the x^ function around the minimum. That is, 

F 

X'(P) - x'(Po) ^ Ax'(x) = x^>l X = em = Y. ^^^^ ' (4) 

where 



i=l 



Xi 



(P - Po)i 
(Po)» 



(5) 



are scaled dimensionless variables, A4 = ADA'" , and ^ = ^-'"x are dimensionless variables in a rotated basis. Here 
A is the orthogonal matrix whose columns are composed of the normalized eigenvectors and T) = diag(Ai, . . . , \f) 
is the diagonal matrix of eigenvalues. Although the full covariance analysis may be carried out without the need to 
diagonalize the matrix of second derivatives M., doing so is both simple and illuminating. For example, the small 
oscillations around the x^-minimum may be represented as a collection of F uncoupled harmonic oscillators. Doing 
so readily identifies the "stiff" and "soft" modes in parameter space. In particular, the soft modes indicate the 
linear combination of parameters that are poorly constrained by the choice of observables included in the quality 
measure. In turn, this suggets the kind of additional physical observables that are required to further constrain the 
model. In the particular case of FSUGold, the softest direction is dominated by the "in-phase motion" of the isovector 
coupling constants gp and Av [7j. This particular linear combination remains largely unconstrained because of our 
poor knowledge of the density dependence of the symmetry energy. 

To obtain meaningful theoretical uncertainties as well as to assess the degree of correlation between two observables, 
one must compute the statistical covariance of two observables A and B. Assuming that the model parameters are 

distributed according to quality measure as exp ( —x^Ai x/2 ) , the covariance of A and B may be written as follows: 



.X ^N ^ dA,,-^ ,, OB J^dA^ .dB 

i,j — 1 "* i—1 

Note that the variance cr^(A) of a given observable A is simply given by cr^(^) = cov{A,A). Also note that, all 
other things being equal, the cov(^, B) is dominated by the softest direction. Finally, the Pearson product-moment 



correlation coefficient — or correlation coefficient for simplicity — is defined as: 



y'va,T(A)va,r{B) 
Further details on the covariance analysis may be found in Refs. [71 IT71E5] . 

III. RESULTS 

The main goal of this contribution is to use a covariance analysis to reliably assess the correlation between the 
neutron-skin thickness of ^"^Pb (i?skin) and a large number of observables of relevance to the structure, dynamics, and 
composition of neutron stars. Moreover, we rely on a covariance analysis to attach meaningful uncertainty estimates 
to our theoretical predictions. In particular, the FSUGold model predicts the following neutron-skin thickness of 
208p|-, vifith properly computed theoretical errors [^i 

i?,i,i„ = (0.2069 ± 0.0366) fm [17.698%] . (8) 

Note that this prediction fits comfortably within the recently determined value by the PREX collaboration [see 
Eq. (fl])]. However, there is a significant difference in the central value that may prove very interesting if the follow-up 
PREX measurement can significantly reduce the error bars without dramatically affecting the central value. 

Given the large number of observables that will be discussed, the present section has been divided into several 
subsections. In particular, each subsection motivates the connection of a given neutron-star observable to i?skin 
and discusses its relevance in constraining the density dependence of the symmetry energy in regions that may be 
inaccessible to laboratory experiments. 

A. Pure Neutron Matter 

Assuming the validity of General Relativity, the equation of state of cold, catalyzed neutron-rich matter is the sole 
ingredient that determines the structure of spherical neutron stars in hydrostatic equilibrium. As such, one of the most 
stringent constraints on the structure of neutron stars comes from the EOS of pure neutron matter (PNM). Although 
PNM remains a theoretical construct, enormous progress has been made in constraining its equation of state at low 
densities. Interestingly enough, most of the progress in this area has been driven by remarkable advances in cold-atom 
experiments that make possible to study the universal behavior of resonant Fermi gases at the unitary limit of infinite 
scattering length. Even though the appreciable effective range of the neutron-neutron interaction invalidates some 
of the powerful arguments associated with resonant Fermi gases, model independent results have been obtained for 
dilute neutron matter [30] . Moreover, the full power of quantum Monte Carlo methods has been used to extend the 
calculations to higher densities [311 132]- By doing so, the calculated EOS — which matches smoothly to the analytic 
results — provides important constraints on nuclear density functionals. Indeed, the impact of such constraints can be 
clearly assessed in Fig. [l] where the equation of state of pure neutron matter for a variety of microscopic approaches 
(see Refs. [31] |31] and references therein) is displayed alongside the predictions from three relativistic functionals. 
Note that although not included in the calibration procedure, the EOS predicted by the FSUGold interaction (solid 
blue line) appears consistent with most of the microscopic approaches. Also note that the required "softening" of 
the EOS relative to the more traditional RMF models (such as NL3) was generated by constraining the calibration 
procedure by the dynamics of nuclear collective modes |18j . 

The wide (blue) band in Fig. IT] represents the la uncertainty in the predictions of the FSUGold model. That the 
Icr-band is so wide reflects our poor understanding of the density dependence of the symmetry energy (recall that a 
20% uncertainty was assumed for the slope of the symmetry energy L). In order to reduce such a large uncertainty 
one must constrain L — or equivalently the pressure of PNM at saturation — through an accurate measurement of the 
neutron radius in ^°*Pb. Although the correlation between the neutron radius of ^''^Pb and the slope of the symmetry 
energy L is by now very well established, only recently has a statistically meaningful correlation coefficient between 
such quantities been computed [Tj [17] . We have now extended such a covariance analysis to explore the correlation 
between the neutron-skin thickness of ^°^Pb and the pressure of PNM at half, one, and two times nuclear-matter 
saturation density (see Table |l] and Fig.[2|. 

Given that the symmetry energy is to a very good approximation equal to the difference in energy between PNM 
and symmetric nuclear matter, the slope of the symmetry energy L is closely related to the pressure of PNM at 
saturation density, i.e., P{pq) ~ poL/3. Thus, perhaps not surprisingly, we find correlation coefficients of almost one 
between i?skin and both L and P{po)- Note, however, that the correlation remains as strong for the pressure of PNM 
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FIG. 1: (Color online) Theoretical uncertainties in the energy per neutron as a function of the Fermi momentum for pure 
neutron matter. The shaded blue area are the la theoretical uncertainty associated to the FSUGold model. 



A 


(A) ± AA 


p(y4,i?ekin) 


L(MeV) 


(60.5152± 12.1011) [19.997%] 


0.9952 


Pipo) 


(3.1842 ±0.6349) [19.940%] 


0.9882 


P{Po/2) 


(0.4874 ±0.1721) [35.304%] 


0.9861 


PC^Po) 


(21.8569 ±1.2735) [5.827%] 


0.8016 



TABLE I: Theoretical errors associated to the predictions of the FSUGold model for the slope of the symmetry energy and the 
pressure of pure neutron matter (in units of MeVfm""^) at three different densities. 

at half saturation density P{pq/2) — but deteriorates significantly for the pressure at twice saturation density P{2pq). 
This last fact reflects the insensitivity of flnite-nuclei observables to the high-density component of the equation of 
state. 



B. Neutron-Star Radii 



The structure of neutron stars — particularly the mass-radius relation — depends critically on the equation of state 
of neutron-rich matter. In particular, neutron-star radii provide stringent constraints on the density dependence of 
the symmetry energy. Ideally, one would measure neutron-star radii over a broad range of masses. Indeed, measuring 
neutron-star radii R{M) for a large range of neutron star masses M would allow one to directly deduce the EOS [33] . 
Unfortunately, whereas various neutron-star masses are accurately known |34j , a precise determination of their radii 
does not yet exist. Moreover, constraining the EOS at low densities from the radii of low-mass neutron stars may be 
difficult as these may be very rare. However, the low-density EOS may be constrained from the neutron radii of heavy 
nuclei as these contain similar information [21 I14j . This is because the same pressure that is responsible for supporting 
a (low-mass) neutron-star against gravitational collapse is also responsible for the development of a neutron-rich skin 
in heavy nuclei. Indeed, theoretical predictions suggest a strong correlation between these two observables: the larger 
the neutron skin of a heavy nucleus, the larger the stellar radius. 

To quantify the theoretical uncertainties in neutron-star radii and to assess their correlation to i?skin, we display 
in Fig. [3] the Mass-vs-Radius relation predicted by the FSUGold model. Although the determination of neutron-star 
radii from observations of the luminosity and temperature is both challenging and hindered by uncertainties in models 
of the stellar atmosphere (see Ref. [37] and references therein) significant advances in X-ray astronomy have allowed 
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FIG. 2: (Color online) The correlation coefficient between the neutron-skin thickness of ^"^Pb and the slope of the symmetry 
energy L and the pressure of pure neutron matter at three values of the density. 
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FIG. 3: (Color online) Theoretical predictions for the mass-versus-radius relation of neutron stars calculated within Icr uncer- 
tainty in the FSUGold model. The observational data that suggest very small stellar radii represent Ict confidence contours for 
the three neutron stars reported in Ref. 35.. The two shaded areas that suggest larger radii are la and 2cr contours extracted 
from the analysis of Ref. j36j . 



the simultaneous determination of masses and radii from a systematic study of several X-ray bursters |35[ I36j . Results 
from such studies are displayed in Fig. |3] alongside the FSUGold predictions. We should note that while we believe 
that studies of X-ray bursters will eventually become instrumental in constraining the dense matter equation of state, 
at present they suffer from systematic uncertainties |Mj. For example, the analysis from Ref. [3S] (occupying the 
upper left-hand corner of the figure) suggests very small radii that are difficult to reconcile with the predictions from 
relativistic mean-field models [38 . In contrast, the more recent analysis by Steiner, Lattimer, and Brown ^36; suggests 
significantly larger neutron-star radii that appear consistent with the predictions of the model. Note that the wide 
(blue) band in Fig. Is] denotes our theoretical uncertainties in the predictions of the stellar radii. We observe that 
whereas the stellar radius is fairly well constrained for large-mass neutron stars, radii of low-mass stars display large 



theoretical errors. For example, the radius of a M = IAMq neutron star is determined with a 3.6% uncertainty 
(Tableln]). Yet the error grows by a factor of three for a M = 0.6Mq neutron star. Such behavior is dominated by 
the large uncertainty in L. Whereas the radius of a heavy neutron star is sensitive to the equation of state at both 
intermediate and high densities, the radii of low-mass neutron stars is sensitive to the range of densities probed in 
heavy nuclei. Indeed, FSUGold predicts a central density in a M — 0.5Mq neutron star of 0.21 fm~ , only slightly 
higher than saturation density. Due to the large sensitivity of neutron stars radii to L, a correlation between the 
neutron-skin thickness of ^°^Pb and stellar radii has been established [H [31 [T3]. Moreover, we expect a stronger 
correlation between i?skin and the radii of low-mass neutron stars than between i?skin and the radii of massive stars. 
This fact is clearly borne out by the results displayed in Table[TTl and Fig.|4J 
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(A) ± AA 


p{A, 7?skin) 


Ro.6 


(13.9785 ± 1.5183) [10.862%] 


0.9953 


Ro.s 


(13.5204 ± 1.0446) [7.726%] 


0.9931 


Ri.o 


(13.2439 ± 0.7776) [5.872%] 


0.9866 


Rl.2 


(12.9864 ±0.5964) [4.593%] 


0.9770 


RlA 


(12.6568 ± 0.4603) [3.637%] 


0.9486 


Rl.G 


(12.1038 ±0.3881) [3.206%] 


0.8361 



TABLE II: Theoretical errors associated to the predictions of the FSUGold model for the radii (in km) of neutron stars of 
various masses. Note that the subscript indicates the neutron-star mass in solar masses. 
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FIG. 4: (Golor online) The correlation coefBcient between the neutron-skin thickness of ^°*Pb and the radius of neutron stars 
of a variety masses. Note that the subscript indicates the neutron-star mass in solar masses. 



Direct Urea Proeess 



Neutron stars are born very hot (with a temperature of about 10^^ K) and cool rapidly via neutrino emission 
through the direct "Urea" process involving neutron beta decay and electron capture 



n 



p + e 



n + Ve 



(9a) 
(9b) 



After most of the "neutronization" process is complete, the standard cooling scenario assumes that neutrino emission 
proceeds through the modiGed Urea process: 



n —^ n + p + e 



(10) 



However, given that the modified Urea process requires the presence of a bystander nucleon to conserve momentum, 
this process may be millions of times slower relative to the direct Urea rate [41]. The putative transition into the 
significantly slower modified Urea phase is solely based on the assumption that the proton fraction in the stellar core 
is below the ^15% required to conserve momentum at the Fermi surface. However, such an assumption is suspect as 
the proton fraction is determined entirely by the density dependence of the symmetry energy, which remains poorly 
constrained. In particular, a stiff symmetry energy, namely, one that increases rapidly with density, favors large proton 
fractions and this may facilitate some enhanced cooling. Note, however, that unlike other enhanced-cooling scenarios 
that may involve the presence of exotic particles in the core, this enhanced-mechanism is not exotic as it only involves 
standard particles and a relatively stiff nuclear symmetry energy. Also note that recent X-ray observations by the 
Chandra observatory seem to suggest that some neutron stars may indeed require some form of enhanced cooling 
(see [12] and references therein). 

Given that the proton fraction in neutron-star matter is controlled by the density dependence of the symmetry 
energy, a strong correlation between the onset of the direct Urea process and the neutron-skin thickness of ^"^Pb has 
been established [15] . Momentum conservation at the Fermi surface demands that the sum of the Fermi momenta of 
the protons and electrons must be greater than or equal than the neutron Fermi momentum. That is. 



rvp ^ /Cp 



(11) 



The onset for the direct Urea process {kp = kp + kp) together with the condition of charge neutrality is sufficient to 
determine the threshold proton fraction: 



T^Urca 
P 



l + {l + xl/^f 



.1.0.11 
=1 9 



(12) 



In the above expression Xc is the electron-to-proton fraction and the arrow indicates the model-independent limit of 



1/9 for the case of a vanishing muon fraction (x^ = 1- 



= 0). For the realistic case of a non-zero muon fraction, the 



threshold proton fraction is model dependent and increases slightly to about yUrca < q ^5. In particular, FSUGold 
predicts a threshold proton fraction of yUrca _ 0.137. That is, the FSUGold model predicts that if the proton 
fraction within the neutron star exceeds this threshold value, then enhanced cooling via the direct Urea process is 
possible. Note that this threshold proton fraction is reached at a density of about three times saturation density which 
corresponds to the central density of a Murca = l-3OM0 neutron star (see Table HI). This suggests — in the particular 
case of the FSUGold model — that any neutron star with a mass below 1.30 Af© that displays enhanced cooling is 
likely to contain an exotic core. For comparison, we should mention that for the significantly stiffer NL3 equation of 
state, the threshold proton fraction (11^'''^'* = 0.129) is reached at the significantly lower density of purca = 0.205 fm~ . 
This would imply that any neutron star with a mass in excess of 0.84 Mq will cool rapidly by the direct Urea process. 
Finally, we note that regardless of the model, the muon-to-proton fraction at the threshold density is fairly significant; 
of the order of 30-40%. 
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(A) ± AA 


p{A,Rskin) 


PUrca 


(0.4668 ±0.1324) [28.359%] 


-0.9928 


Murca /Mq 


(1.3012 ± 0.2658) [20.427%] 


-0.9927 


-i^Urca 


(0.1367 ±0.0019) [1.421%] 


-0.9927 


YA^Po) 


(0.1064 ±0.0138) [13.000%] 


±0.9906 


Ypipo) 


(0.0609 ± 0.0055) [9.055%] 


±0.9166 


YApom 


(0.0346 ±0.0051) [14.651%] 


-0.9063 



TABLE III: Theoretical errors associated to the predictions of the FSUGold model for various neutron-star properties of 
relevance to the direct Urea process. The threshold density is given in units of fm""^. 



To quantify the strong correlation between var ious neutron-star properties at the Urea threshold and the neutron- 
skin thickness of ^"^^Pb, we have listed in Table HI and plotted in Fig.js] the corresponding theoretical errors and 
correlation coefficients extracted from our covariance analysis. As suggested above, we find a strong direct correlation 
between the proton fraction at high densities and i?skin- Note, however, that at low densities the proton fraction and 
Rskin become aiiti-correlated, as a stiff symmetry energy goes to zero faster than a soft one. Moreover, as expected, the 
Urea threshold density, stellar mass, and proton fraction display a strong (inverse) correlation to i?skin- In particular. 
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FIG. 5; (Color online) The correlation coefficient between the neutron-skin thickness of ^"^Pb and various neutron-star prop- 
erties of relevance to the direct Urea process. 

we note that the strong anti-correlation between yUrca g^^^ i?skin emerges entirely because of the presence of muons 
in the star. Indeed, without muons Y^'^'^^ is fixed at 1/9 so all of its derivatives with respect to the model parameters 
will vanish and, thus, so will the correlation coefficient [see Eq. Q]. Finally, we have shown that a thin neutron skin 
in ^"^Pb, i.e., a very soft symmetry energy, requires a large neutron-star mass for the onset of the direct Urea process. 
Thus, we believe that the observation of enhanced cooling in low-mass neutron stars (M < Afurca) provides one of the 
most promising indicators of exotic states of matter residing in the stellar cores. 



D. Core-Crust Transition 



Neutron stars have a solid crust above a uniform liquid core. The structure and composition of the solid crust 
remains a source of significant debate and of considerable interest [43H54] . Moreover, the crust is believed to be of 
critical importance in a variety of fascinating astrophysical phenomena, such as pulsar glitches, giant magnetar flares, 
and the emission of gravitational waves [55H58J . The phase transition from the solid crust to the uniform liquid core 
depends on the properties of neutron-rich matter. As neutron-rich matter becomes dilute, the uniform ground state 
becomes unstable against small-amplitude density fluctuations O [14] . That is, it becomes energetically favorable for 
the system to fragment into high-density clusters (i.e., nuclei) embedded in a dilute neutron-rich vapor. Yet details 
of the transition depend sensitively on the proton fraction — and thus on the density dependence of the symmetry 
energy — at low densities. Whereas the proton fraction at high density displays a direct correlation to the neutron- 
skin thickness of ^°®Pb, the proton fraction at low densities is anti-correlated to i?skin- For example, the correlation 
coefficient between Yp{pQ/2) and i?skin is both large and negative (see Table III). Given that the symmetry energy 



represents the energy cost in departing from the isospin symmetric limit, a stiff symmetry energy falls rapidly to 
zero at low densities and, hence, tolerates a larger isospin asymmetry (i.e., smaller Yp) than their softer counterparts. 
Thus, a lower proton fraction typically implies a low transition density from the solid crust to the liquid core. This 
suggests an inverse correlation: the thicker the neutron skin thickness of ^'^^Pb, the lower the core-crust transition 
density [Z\ 

Several approaches — of both a microscopic and thermodynamic nature — have been used to determine the insta- 
bility of the uniform ground-state against cluster formation [3l [Ml [SH [59H6T] . In the present study we rely on the 
RPA stability analysis described in Ref. '14' to compute the baryon density, proton fraction, pressure, and energy 
density at the core-crust interface. Results from the covariance analysis are listed in TablelTV] and displayed in 
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A 


(A) ± AA 


p{A, i?skin) 


Pt 


(0.4020 ±0.1071) [26.640%] 

(0.0351 ±0.0069) [19.711%] 

(71.5337 ±5.3747) [7.514%] 

(0.0755 ± 0.0056) [7.369%] 


±0.9474 
-0.9260 
-0.9207 
-0.9203 



TABLE IV: Theoretical errors associated to the predictions of the FSUGold model for various neutron-star properties of 
relevance to the transition between the solid crust and the uniform liquid core. The transition pressure and energy density are 
given in units of MeVfm"'^ and the transition baryon density in fm~^. 
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FIG. 6: (Color online) The correlation coefficient between the neutron-skin thickness of ^"^Pb and various neutron-star prop- 
erties of relevance to the transition between the solid crust and the uniform liquid core. 

Fig. [6j For comparison, the stifFer NL3 equation of state — which has a significantly smaller symmetry energy at 
sub-saturation densities than FSUGold — predicts the following significantly lower values: pt = 0.052 fm~ , 11* =0.015, 
Pt = 0.212 MeVfm"^, and £t = 48.960 Me Vfm~^. As suggested earlier, models with a stiffer symmetry energy and 
thus thicker neutron skins, display a "delayed" transition from the uniform core to the solid crust. Thus, we find a 
strong anti-correlation between i?skin and the various stellar properties at the core-crust interface [3]. Yet, the tran- 
sition pressure P^ behaves in an interesting an unique way. First, in contrast to other observables, Pt appears to be 
directly correlated to i?skin with a large and positive correlation coefficient of about 0.95. Second, the mere existence 
of a strong correlation appears to be at odds with earlier studies that suggest a weak correlation between Pt and 
-Rskin [El ISH IM] . We attribute the present intriguing result to the fact that within the realm of a covariance analysis, 
the correlation coefficient between two observables is obtained by generating model parameters that are distributed 
according to the quality measure [see Eq. Q]. That is, models that differ significantly from the optimal (FSUGold) 
model carry little weight. In contrast, the weak correlation between Pt and Pskin suggested in Ref. [IS] was obtained 
through a systematic (perhaps even ad-hoc) variation around the optimal model. Although the covariance analysis 
implemented here provides the proper statistical measure of correlation between two observables [29] . a covariance 
analysis can not assess systematic errors associated with the limitations of a given model. Thus, whereas the FSUGold 
model — with its soft symmetry energy — predicts a strong positive correlation between Pt and Pskin, a model with a 
stiffer EOS may predict exactly the opposite (see Fig. 6 of Ref. [16j). 



E. Stellar Moment of Inertia 



Although the general-relativistic expression for the moment of inertia of a neutron star is fairly complex — even in 
the so-called slow-rotation approximation [621 [S3] — on simple dimensional grounds it must scale as the product of 
the stellar mass times the square of its radius. Thus, the possibility of measuring the moment of inertia to even a 
10% accuracy may provide stringent constraints on the nuclear equations of state [64H67] . The prospects for such a 
measurement improved significantly with the discovery of the binary pulsar PSR J0737-3039 [SHI ISS]- The ffist ever 
discovered double pulsar, PSR J0737-3039 exhibits characteristics that have enabled the accurate determination of 
several pulsar properties (such as the orbital period of the binary, both pulsar masses, and both spin periods) and 
have resulted in some of the most precise tests of Einstein's general theory of relativity. 
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Motivated by these facts, we have recently studied the sensitivity of the stellar moment of inertia to the equation 
of state |l6j . However, we only found a mild sensitivity of the total moment of inertia to the underlying EOS. Thus, 
we redirected the main focus of such contribution to the crustal component of the moment of inertia. Several reasons 
prompted this choice. First, constraints on the EOS may be imposed from an analysis of pulsar glitches in the Vela 
pulsar that place at least a 1.4% of the total moment of inertia in the solid crust 55_ . Second, the crust is thin and the 
density within it is low, so fairly accurate analytic expressions for the crustal moment of inertia have been developed 
in terms of two stellar properties that are highly sensitive to the EOS: (a) the radius of the uniform core and (b) the 
transition pressure at the core-crust interface 16, 55j [56] . Finally, given the strong correlation just found between 
the core-crust transition pressure Pt and i?skin, one expects the emergence of a similar strong correlation between the 
crustal moment of inertia and i?skin- Note, however, that our expectation of a strong correlation is solely based on 
the results of the covariance analysis. If instead one relies on a systematic variations around an optimal model, no 
strong correlation was found |16| . 



A 


(A) ± AA 


p{A, Rskin) 


7-cr 

rcr 

-'1.4 

^0.8 

/l.4 


(8.7777 ±2.5612) [29.178%] 

(5.8988 ± 1.4055) [23.827%] 

(7.4067 ±0.3204) [4.326%] 

(14.7660 ± 0.3437) [2.327%] 


0.9781 
0.9619 
0.9299 
0.5192 



TABLE V: Theoretical errors associated to the predictions of the FSUGold model for the crustal (in units of 10' 
total (in units of IC*^ gcm^) moment of inertia of a 0.8 and 1.4 solar-mass neutron star. 
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FIG. 7; (Color online) The correlation coefficient between the neutron-skin thickness of ^"^Pb and various neutron-star prop- 
erties of relevance to the crustal and total moment of inertia. 



In Table[V|we list the crustal and total moment of inertia for neutron stars of 0.8 and 1.4 solar masses; the correlation 
coefficients are also displayed in graphical form in Fig. [7] As expected, the crustal moment of inertia — being sensitive 
to the core-crust transition pressure Pt — displays a strong correlation to i?skin- Note that the large error attached to 
this observable reflects the large theoretical uncertainty associated with i?skin — and ultimately with L. Given that the 
central density for a low-mass neutron star is only slightly larger than saturation density, the correlation between the 
total moment of inertia of a 0.8 Mq neutron star and i?skin remains strong. However, we find a significant deterioration 
in the correlation between i?skin and the moment of inertia of a 1.4 solar-mass neutron star. 

We close this section by collecting in Fig. [8] all correlation coefficients displayed earlier between i?skin and the large 
number of neutron-star observables. Moreover, we display in a color-coded format in Fig.|9] correlation coefficients 
for 16 observables (i.e., 120 independent pairs) of relevance to the structure, dynamics, and composition of neutron 
stars. As shown earlier, most of these stellar observables display a strong correlation — or anti-correlation — to the 
neutron skin thickness of ^°®Pb (first column/row in Fig.^. As mentioned above, a notable exception is the total 
moment of inertia of a canonical 1.4 Mq neutron star. Of course, not every neutron-star observable is sensitive to the 
density dependence of the symmetry energy. The maximum neutron-star mass Mmax? with a correlation coefficient 
of only p(M,„ax, -Rskin) = 0.0163, provides a particularly clear example. Note, however, that M^ax yields one the best 
constraints — if not the best — on the high density component of the equation of state. 
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FIG. 8: (Color online) The correlation coefficient between the neutron-skin thickness of ■^°*Pb and the large number of neutron- 
star properties discussed in the text. 



IV. CONCLUSIONS 



The successfully commissioned Lead Radius Experiment (PREX) at the Jefferson Laboratory has provided the first 
largely model-independent determination of neutron-skin thickness of ^"^Pb [12]: i?skin = Rn — Rp = CSSl^Ql^g fm. 
Given that the determination of the neutron radius of a heavy nucleus is known to have strong implications on the 
structure, dynamics, and composition of neutron stars pi-Hl [Tll - [roil70| . a detailed covariance analysis was implemented 
to quantify the correlations between -Rskin and a variety of neutron-star observables. Moreover, meaningful theoretical 
uncertainties were provided for all predicted observables. We note that although the significant impact of a measure- 
ment of i?skin on neutron-star properties has been known for almost a decade, our work quantifies for the first time 
these correlations on the basis of a detailed covariance analysis. We stress that the covariance analysis employed here 
represents the least biased and most comprehensive tool to uncover correlations between physical observables [71 117) . 

In the present study we have computed correlation coefficients between i?skin and the following quantities of direct 
relevance to the physics of neutron stars: (a) the equation of state of pure neutron matter, (b) stellar radii, (c) the 
onset of the direct Urea proces, (d) the core-crust phase transition, and (e) the total and crustal moment of inertia. In 
addition to their intrinsic importance, these observables are interesting as they are sensitive to the density dependence 
of the symmetry energy over a wide range of densities. For example, whereas the onset of the direct Urea process is 
sensitive to the symmetry energy at high densities, the core-crust transition probes the symmetry energy at densities 
of about a third to a half of nuclear-matter saturation density. 

Although most of the correlations follow the expected trend, two observables — both strongly correlated to -Rskin — 
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FIG. 9: (Color online) Color-coded plot of the 120 independent correlation coefficients between 16 physical observables of 
relevance to the structure and dynamics of neutron stars. 



deserve a special mention. These are the threshold proton fraction for the direct Urea process Y^'^'^^ and the core-crust 
transition pressure Pt. First, the strong anti-correlation (of about —0.99) between yUrca ^^^ i?sijin came as a surprise 
because the sensitivity of yUrca ^^ ^^^ density dependence of the symmetry energy is due entirely to the presence 
of muons in the stellar core. Indeed, in the absence of muons yU'''^^ is fixed at 1/9 and the correlation coefficient 
vanishes. Second, the strong direct correlation between Pt and -Rskin (of about -f 0.95) is surprising because no such 
correlation was observed in previous studies }16, :59i, :61s . Note, however, that those earlier studies do not rely on a 
covariance analysis, but rather explore possible correlations by using either a large number of nuclear models [SSI ISI] 
or models that are systematically varied around an optimal one |16] . Although the covariance analysis implemented 
here provides the proper statistical measure of correlation between two observables |29j , a covariance analysis can not 
assess systematic errors associated with the limitations of a given model. Thus, the implementation of a covariance 
analysis using other nuclear functionals is both highly desirable and strongly encouraged. 

In connection to the follow-up PREX measurement that aims to achieve the original 1% precision goal, we regard it 
of a great value, as few accurately calibrated models — if any — predict such a large (central) value for -Rskin- Moreover, 
such a large neutron skin — which would imply a fairly stiff symmetry energy — will have widespread repercussions in 
the physics of neutron stars. It is worth mentioning, however, that a 10% uncertainty in the determination of the slope 
of the symmetry energy appears to require a more stringent measurement (at the 0.5% level) of the neutron radius of 
^•^^Pb [HIITT]. Ultimately, a determination of the density dependence of the symmetry energy is likely to require a 
multi-pronged approach [55]. First, from the theoretical perspective, we expect that powerful arguments based on the 
universality of dilute Fermi gases in the unitary regime will continue to shed valuable insights into the behavior of pure 
neutron matter. Second, in terms of laboratory observables, the dipole polarizability of ^°^Pb — measured recently 
with unprecedented accuracy |72j — provides a unique constraint on the density dependence of the symmetry energy 
and an excellent complement to the neutron-skin thickness of ^"^Pb ^3J. Finally, enormous advances in both land- 
and spaced-based observatories have started to impose significant constraints on the equation of state. In this work we 
have demonstrated the existence of many neutron-star observables that — by virtue of their strong correlation — may 
be used as a proxy for i?skin- A particularly interesting alternative is the radius of a low-mass neutron star, as its 
central density is close to nuclear-matter saturation density. Indeed, the radius of 0.8 Mq neutron star, Rq.s, displays 
a strong correlation (of ^0.99) to i?skin- Thus, a 5% measurement of i?o.8 would translate into a ~0.4% constraint 
on the neutron radius of ^''^Pb. Unfortunately, the measurement of stellar radii remains a serious challenge that is 
further complicated by the scarcity of low-mass neutron stars. Given that most well measured neutron-star masses 



14 

lie within a narrow range centered around IAMq, low-mass neutron stars may be difficult to form and therefore may 
not even exist. In such a case, a highly accurate measurement of -Rgkin remains the sole alternative. 
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